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We study the basic thermodynamic and electromagnetic properties of the super- 
conductor described by the negative-?/ Hubbard model ( gap parameter A, critical 
temperature T c , London penetration depth A, thermodynamic critical field H c and 
Ginzburg-Landau correlation length £g-l)- 

The calculations are performed for square (SQ, d=2), simple cubic (SC, d=3) 
and face-centered cubic (FCC, d=3) lattices. We analyze the results as a function 
of electron density n and interaction U using the Hartree-Fock approximation in the 
weak-to-intermediate U case combined with the conclusions from the perturbation 
theory valid in the strong coupling limit. From the equality: A/£ G _ L = l/v / 2 we find 
the boundaries between local and nonlocal electronic behaviour of the U < Hub- 
bard model in the U — n parameter space. Using the calculated values of T c and 
A we compare the 'universal plot' of T c /T™ vs. [A(0)/A m (0)] 2 ( where: T™ denotes 
maximum T c as a function of n and A m (0) corresponds to T™) with the recent exper- 
imental data for various families of nonconventional superconductors and conclude 
that the best agreement can be obtained for intermediate values of the local attrac- 
tion. 

Keywords: negative-?/ Hubbard model, London penetration depth, high-T c super- 
conductivity 
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I. INTRODUCTION 



Among the recently studied superconducting compounds there are several groups of 
materials (copper oxides, fullerides, Chevrel phases, barium bismuthates ...) which share 
some nonconventional features like high T c with relatively low carrier concentration, short 
coherence length together with extremely small £g-l/A ratio and universal dependence of 
T c /T™ vs. [A(0)/A m (0)] 2 Above features suggest that a short range almost nonretarded 
attraction is responsible for pairing in these systems || . One of the simplest effective models 
to describe this situation is the negative-?/ Hubbard hamiltonian. This model constitutes a 
common basis for the description of the superconductors with weak local electron pairing, 
being in many ways similar to the classical BCS systems, and superconductors with strong 
attraction where the local pairs conform to weakly charged hard core bosons @. It has been 
considered as an effective model of superconductivity in the family of cuprates , the 

barium bismuthates |7||1, the fullerides and the Chevrel phases 0. 



In our paper we present some new results of an analysis of the several basic thermody- 
namic and electromagnetic properties of this model for various types of lattice structure. 
We have concentrated on the case of (i) square (SQ, d=2), (ii) simple cubic (SC, d=3) and 
(hi) face-centered cubic (FCC, d—3) lattices, i.e. structures representing (i) the family of 
cuprates, (ii) the Chevrel phases and barium bismuthates and (iii) the fullerides, respec- 
tively. For comparison we also present the results obtained for a model rectangular density 
of states T>(e). To calculate quantities of interest we use the (broken symmetry) Hartree- 
Fock approximation (HFA) which is known to give credible results at T = K as concerns 
energy of the ground state, the energy gap, the chemical potential as well as the collective 
excitations in the whole interaction range, interpolating smoothly between the weak and 
strong-?/ limit, where it matches the results of the perturbation theory developed from the 



zero-bandwidth limit [14]. On the other hand the HFA leads to qualitatively erroneous 



results for the the thermodynamic critical field H c and Ginzburg-Landau correlation length 



£ G _ L in the strong attraction limit. This failure is due to the fact that the HFA greatly 
overestimates the energy of the normal state which is used in a standard calculation of H c 
and next £ G _ L . In order to discuss the behaviour of H c , £ G _ L and Ginzburg ratio k = | in the 
whole interaction range we have completed the results of the HFA calculation with the ones 
of the perturbation theory. We will also point out a possible relevance of our results to the 
recent experimental data for the short-coherence length superconductors. 



II. GENERAL FORMULATION 

We consider the negative-?/ Hubbard Hamiltonian in the hypercubic lattice 

H = tJ2 (e^^ct a c^ a + H.c.) - \U\ £ n^n Al , (1) 

fhSo m 

where t is the hopping integral, 5 is a unit lattice vector and the exponential (Peierls) factors 
in (JJ) account for the coupling of electrons to the magnetic field via its vector potential A(f). 
$ mm+< 5 = — 4-A(m)5, and e is the electron charge. 

For weak vector potential the expectation value of the Fourier transform of the total 
current operator in a direction (a = x, y, z) can be obtained from the linear response theory 
T|,0 as 



J a (q,u;) = N^^[5 a pK di « + K p J 1 (q,u J )]Ap(q,u;) . (2) 

Here we have separated the total response kernel K a /3(q,uj) into a diamagnetic part and a 
paramagnetic part. The diamagnetic contribution evaluated within HFA is 

K am = 9 - — — > <cjct >cosk Q = 
h 2 c 2 a ± N 

ka 

8vre 2 |t| 1 



h 2 c 2 a ± N 

k 



cos k a , (3) 



where Ag = eg — fi, eg = t Y.g e an d ^% = + ^ 2 * s ^ e q. uas ip ar ticle energy with the 
superconducting gap parameter 



A = \U\x = lJ ff-Ei < 44 >= "tv ^fc < 4t c ^fci > The chemical potential ia fi = fl-\U\n/2, 



M 



fc? -fcj. 



iV and n = Jjf denote the number of lattice sites and the number of electrons per site, 
respectively and < n < 2. The values of A and \x are obtained as a solution of the system 



of the selfconsistent equations |L5 



A = ^Etanh(^/2)/^ , 



n- 



^E[!-^tanh(^/2) / E k 



(4) 
(5) 



whereas the free energy per site of the superconducting (SS) phase and the normal (NO) 

phase is given by: 

_pss 



N 



=Mn-l) + g-Mn»-^E^ooah^ 
= // NO (n - 1) - ^n 2 - — E In 2 cosh 



AT 



(6) 
(7) 



where /i NO is determined from: 



1 



n - 



- E 1 - tanh 



0A» 



iV 



(8) 



The expression for the paramagnetic part of the kernel reads 

1 167re 2 t 2 
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+ 
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IE 



(9) 



V w " (-Effi-y + + <e a; + £fc + i%_ - 
In the London superconductors the magnetic field penetration depth A(T) is determined 
by the sum of diamagnetic and paramagnetic part of the total kernel in the static limit 
A(T) = lim gH _ [-^* a - Kl™(q x = 0, q y , q z = 0,u = O)]" 1 / 2 0JT7I . 
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At T = K the paramagnetic part of the kernel becomes important in determining A 
when the correlation length becomes greater than the penetration depth and we deal in 
this case with nonlocal (Pippard) superconductor. This situation is common in many low- 
To systems. The short -coherence length materials, including the high T c superconductors 
represent the opposite, i.e. the London limit. In the latter case the ground state penetration 
depth is determined entirely by the q — > limit of the kernel where the paramagnetic part 
of the kernel vanishes and A is given by: 

A= - 1 . (10) 
^/-K dia 

In our calculation of A we restrict ourselves to the London limit having in mind the 
properties of the systems of interest and finaly determine the area in the U-n parameters 
space where the local approximation may be valid. The value of the penetration depth 
calculated in this way is qualitatively good both in the weak and strong-?/ limits, in the 
latter case approaching the results of the perturbation theory, as will be shown below. 

Using the value of the penetration depth and the difference of the free energy between 
normal and superconducting phase one is able to determine the thermodynamic critical field 
H c and the Ginzburg-Landau correlation length £ G _ L as 

H 2 C (T) F NO (T) - T SS (T) 



8n Na ± a 



2 



(11) 



UL ~2n^XH c ' (12) 
where $ — , an d to obtain the estimations for the critical fields H cl ~ ~ ^ and 

H c2 ~ 2^T7 ' where K = li~l' At T = ' l N 1 = E o = jf( H )(T^)- Within HFA the expressions 
determining Eq ° and Eq S are given by: 



E N ° = ^^ = E« o (0)--\U\n* , (13) 



^o ss =^r 2 =^ K s (o)- M - i l^ 2 , (i4) 

where 



-^NO(SS) ~ e fc( C /L C fc<r)NO(SS) , (15) 

kcr 

is the average value of the kinetic energy term in the NO(SS) phase. 

The HFA calculation of the energy of the ordered state at T = K are reliable for any 
U || and Eq S reduces correctly to the exact value: Un/2 in the zero bandwidth limit where 
the electrons form a system of on-site pairs and singly occupancy of sites is prohibited. 
On the other hand the HFA energy of the normal phase in this limit is equal to Un/4 
instead of Un/2 thereby producing the incorrect energy difference of Un/4. In order to 
remedy this inconsistency, at least in part, we resort to the results of the perturbation 
theory in the calculation of for \U\ ^> t. The hamiltonian ([I]) reduces in this case to the 
pseudospin-1/2 Heisenberg model working in the subspace of states with no singly occupied 
sites 
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h = £ <M + - \ £ h ti*js-** m + H-c) 

+ EM 5 ^-f ZJ ' (16) 

where 2p^ = (n^ + — 1), p/" = c^cj, = 2t 2 /|[/| and Z denotes a number of nearest 
neighbours. The hard-core boson operators: pj , p| satisfy the commutation rules of the 
s=2 operators. The electron number condition is: 

"=^E<2pf+l> , (17) 



For xq = (p + ) = t? J2i(pf) 7^ the hamiltonian fll6| ) describes the superconducting state and 



the MFA result for the Eq S is equal to the corresponding HFA one up to terms of order 
t 2 /\U\ inclusive (see Sec. 3). The MFA calculation of the ground state energy for the normal 

state, a:o = 0, gives: 



■n,,_(H>_ l |rrl „ , 1 



E™ = ^ = --\U\n + -ZJn(n-2) , (18) 

Using this expresion in the analysis one should keep in mind that the MFA neglects the 
intersite correlation of fluctuation of the pseudospin operators what can lead to the error of 
order of J = 2t 2 /\U\. 



III. ANALYTICAL RESULTS FOR WEAK AND STRONG U CASES AT T = K 

The system of the selfconsistent equations ([|,|5]) can be approximately solved in the 
limiting cases of the weak and strong attraction limits and here we present these results: 

A) The weak-?/ case 

Provided that the chemical potential is not located too close to van Hove singularity of 
the density of states (DOS) function T>(e), the gap parameter at T = can be obtained 
from equation (||) as: 
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^V^exp (gL) exp (-^) , ( 19 , 

where: .D+ —D2 — p, D_ = D\ + /2 with —D\,D2 denoting lower and upper boundary 
of the band, respectively; /[C^] is a functional defined by 

I[F}= [ D2 design(e-H)r(e) , (20) 

J —D-y 

is a function of electron density, which depends on details of DOS and is defined 

by 

P(e )-P W 

e-/i 

In Eqs.(^9|),(|2T|) the chemical potential may be approximated by its value in the normal 
state, fl NO and it is determined by: 



1 = - f° 2 GfcP(e)sign(e-/2 N ) . (22) 

The critical temperature in this limit is given by: 

k B T c /A « 1/1.76 , (23) 

which is the well-known weak coupling BCS ratio. 

Analogously, performing the weak coupling expansions of Eqs. (|3], ^ [TTJ) at T = one 
can find the analytical expressions determining A, /2 SS , H c , £ G _ L and k. The London 
penetration depth is calculated as 

' ~K^ = -^-\^) , (24) 



n 



A 2 (0) ~ h 2 c 2 a ± Z 

where 
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£ s K s (0)-£* o (0) 



A 2 

W\ 



A 2 j ^NQ-^ 



X 



X 



(£> 2 -/i) 2 P1+//) 2 



+ 



2(D 2 -/2) 2 



(25) 



where 



#fi(<0 



e — [i 



G (e)= I de'V{e')e' 

-Di 



Hie) 



^o(e)= 



(26) 



e - /i J-Di 

and /2 = /2 no . Obviously, for \U\/B — > £^(0) approaches the value of the band 
energy in the normal state. 

In a case when we can neglect variation of T>(e) and D\ = D2 = D Eq . (|25|) reduces to 

A 2 A 2 



£ S K S (0)=-2>(/iNO 

whereas // ss and H c (0) are given by 



n 2 - n 2 J 



(27) 



/4 = /4o 1 + 



D 2 - n 2 ' 

where /i NO is given by Eq. (|2"2"D , and 

tf 2 (0) Eq° — £q s 1 A 2 X>(/i) 



(28) 



57T 



a i a; 



(29) 



2 a ± a 2 

The weak coupling formulas for £ G _ L and ft = A/£ G _ L follow directly from Eqs.(|12"D, (p4|), 
(^9|). In the case of X> 1 (/x) = and Di = D 2 = D their explicit forms are: 

On / 1 \ 



e 



G-L ■ 



D 2 -m Q , i 



i 



2v/2Z^ A 2 ' 2 |tf|D(£) 2^^ \\U\V(ji) 



(30) 



A 2 D(^ 



7TV7T an 



exp(-l/|t/|X>(/2 NO )) 
(D 2 - m o )V(fi NO ) 



A! 



7TV7T a,. 



(31) 
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B) The strong-?/ case 

In this case the integrands in Eqs.(^J^) can be expanded in series with respect to t/U 
and we have: 



A 1 

x = 



t\ 2 _ / t x " 



wr-^ 2 f-^\u) + 165Ls \w\) ' (32) 



fiss = - l ^(l+4L 2 (t/U) 2 -4L 3 (36-l/6)(t/\U\) 3 ) , (33) 



where 5=1 — n, /x NO is given by (^2|) and the lattice sums L n represent the moments 
of DOS and are given by: 

Ln= T,Js 1+ ..J n ,d ■ (34) 

5i...S n 

One can note that odd moments vanish for the alternating (i.e. SQ and SC) lattices 
but not for the FCC ones. For alternating lattices L 2 = Z, L 3 = and for FCC lattice 
L 2 = 12 and L 3 = 48. The corresponding formula for A obtained by expanding of Eq.@ 
reads 

^ = 4^^[2(l-5 2 )^(t/f/) 2 -65(l-5 2 )L 3 (t/|f/|) 3 ] . (35) 
X z n c z a i Z 



In the structures with asymmetric DOS (like FCC) the odd terms in t/U from Eq . (p5| 
are again nonzero. From (|T3| ) and ( |T4|) we obtain: 



K° = -^E^sign(e-/, NO )-^|n 2 , (36) 



\U\ t 2 t 3 

E s s = - l -^n-n(2-n)L 2 — + 2n(l-n)(2-n)L 3 — , (37) 



Hi ES°-E&* 1 [\U\ , . , ,zt 2 \ ^(t 3 \ 



n 



It is worth while to compare the strong coupling results of HFA with the ones obtained 
for the effective pseudospin model Hamiltonian (0). Recently, the electromagnetic proper- 
ties of that model have been analysed in and below we only quote the expressions 
derived within MFA (the RPA treatment, taking into account quantum corrections, yields 
qualitatively similar results, except the low density limit fl9fl ): 

X - 



\\j<2-n) , (41) 



Ef = - l -\U\n-\j Q n{2-n) , (42) 



2' 1 2 



Hi 1 J n(2 - n) 



In 4 a 3 
2ne 2 



(43) 



A- 2 =^-Jn(2-n) , (44) 
n c z a 



where Jq = ZJ, J=jjj\ an d e = 2e, a ± = a ll =a, 



£g-i=^ , (45) 



(46) 



2Z 



£g-l 2ny/n^JaJn(2 — n) 

and Eq° is given by Eq.([18|). Up to terms of order t 2 /\U\ the HFA expressions for xq 
(Eq.(||)), Ef (Eq.fl37|)) and 1/A 2 (Eq.(]3|)) are equal to those given by Eqs.fln]), (||) and 
(|4]). On the contrary the HFA results for H c , £ G _ L and k are qualitatively erroneous due to 
the great overestimation of Eq° by HFA in this limit (comp. Eq . (|36|) and Eq.(|i~8|)). Thus, 
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to get realistic values of H c , £ G _ L and k for arbitrary \U\/B one should combine the results 
obtained using Eq° given by Eq . (|T3l) (reliable for \U\/B < 1) with those obtained using Eq° 
given by Eq.flilf) (reliable for \U\/B > 1). This will be done in the next section, where we 
will also compare some of the predictions given above with the results of numerical solutions 
ofEqs.®-©, ©-(0). 

IV. RESULTS OF NUMERICAL SOLUTION AND DISCUSSION. 

Figs. la and lb show the numerical plots of T c (n), A(n) and X~ 2 (n) for SC and FCC 
lattices, respectively (the corresponding plots for SQ lattice are given in Fig.l of Ref. PHI), 
whereas Fig.lc presents the concentration dependence of A and X~ 2 (n) for the model rect- 
angular T>(e)\ 

V(s) = l/B for - B/2 < e < B/2 = D , otherwise 0, (47) 

where B is the effective bandwidth. All these plots were made for the values of U from 
the weak coupling regime (\U\/B = 0.1) where the effects of the density of states, T>(e), 
are most clearly seen. In this regime the superconducting critical temperature T c and the 
T=0 gap parameter A increase rapidly with p, approaching van Hove singularities in T>(e), 
in agreement with a modified BCS expression Eq. (|T9|) . This is the reason why T c for SQ 
and SC lattices is symmetric with respect to n—1 and so strongly peaked for SQ lattice at 
n=l while the position of T c peak in the FCC structure moves with increase of \U\ from 
the vicinity of toward n—1. For the FCC lattice the effect of density of states on the 
penetration depth consist mainly in [/-dependent shift of the maximum of 1/A 2 from n = 1 
towards n < 1 (Fig. lb). The inverse squared penetration depth, being proportional to the 
bandwidth for weak-U decreases like ~ Zt 2 /\U\ in the strong coupling limit. The plot of 
this quantity as a function of n for a few values of \U\/B is shown in Fig. 2 (SQ lattice). 
Examples of the evolution of the thermodynamic critical field H c and of the correlation 
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length £ G -L) with \U\ and n are shown in Fig.3,4. Note the steep decrease of £ G _ L with n in 
a relatively narrow range of n values for the FCC lattice and substantial reduction of £ G _ L 
near n = l for the SQ lattice (Fig. 5). These features can be understood on the basis of the 
A(n) behaviour presented in Figs 1,2. 

The substantial variation of £ G _ L with n for \U\/B < 1 is largely due to the strong n 
dependence of H 2 which in the weak coupling limit is proportional to A 2 V(fi) (cf. Eqs.(p9|), 

«■ 

With increasing \U\ the n-dependence of the correlation length is less pronounced and 
in the strong coupling regime, \Zt/U\<^l, £ G _ L calculated with in HFA goes like ~ \Zt/U\ 
being almost independent on n. However, for large \U\ the energy of the normal phase, 
entering (|ITD , ( |I2"D is drastically overestimated in the HFA (comp. Sec. 3). 

To get more realistic values of H c and £ G _ L for \U\/B>1 we have made an estimation of 
£ G _ L and H c exploiting perturbation approach in the calculation of the energy of the normal 
phase (Eq.(^)) and the numerical results for SQ lattice and n — 1 are shown in Figs 3,4. In 
the \U\/B 3> 1 limit the analytical results for H c , £ G _ L and k can be obtained in this way 



by using Eqs.Qlq, |37| , |35| and 



Hi ES°-m* 1 (I . . t 2 . w t 3 \ 

8f = J ^ = AT {- 2 n(2-n)L w -2n { l-n) { 2-n)L 3 - 2 \ , (48) 



,G_L V2Z' 



3(1 -n)L z 4, 



\Z-A{l-n)L^ 



$o ZJEZ 1 J z -4(l -n)L 3Wl 



(49) 



(50) 



27T0F o„ ^(2-n)^ Z- 3(1 -n)L 3 ± 

They coincide (to the second order in t^t) with the MFA results for the pseudospin model 
(0) (cf. Eqs.fl, (H), ©). 

As we see from Figs 3a and 4, for \U\/B > 1 H c decreases with increasing \U\/B and 
for \U\/t ^> 1, H 2 ~ t 2 /\U\, whereas £ G _ L tends to a fixed value a/\ / 2Z. 
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Finally let us conclude the evolution of the Ginzburg ratio k = A/£ G _ L with \U\ and n. 
Concentration dependence of this quantity is plotted in Fig 6. Notice universal behaviours (i) 
K ~ n~ 1 / 2 for n <C 1 (arbitrary \U\/B) and (ii) k ~ [n(2 — n)] -1 / 2 for \U\/t ^> 1 (arbitrary n). 
In the weak-to-intermediate coupling regime (\U\/B < 1) k{ji) is not universal and strongly 
depends on the details of V(e). From the equation A/£ G _ L = l/v2 one can estimate the 
boundaries between local and nonlocal electromagnetic behaviour in the considered model 
fTE|,I7|. In Fig. 7 we plot this boundary in the \U\/B — n parameter space for fixed value 
of t and a. In the low concentration limits (|n — 1| ~ 1) the local behaviour extends up to 
very small values of \U\/B. 

The high-T c systems (not only the cuprates but also the barium bismuthates and the 
fullerides) are extreme type II superconductors with local electromagnetic behaviour and 
with very large Ginzburg ratio estimated to be of order 10 2 |3,21|. Thus, from the equality 
k— 10 2 one can estimate the values of \U\/B and n which could be reliable for these materials. 
Examples of such estimations are given in Fig. 8 and they indicate that best agreement can 
be obtained for intermediate values of the local attraction. 

In order to make some more detailed comparisons with experiment we have evaluated 
the plots T c /T™ vs. (A m /A) 2 , where T™ denotes the maximum critical temperature in the 
0<n<2 range and A m corresponds to T™. For a given \U\fB the data corresponding to SQ 
and SC lattices represent the one-valued function (see Figs.9ab), whereas the ones for the 
FCC lattice (see Fig. 9c) exhibit substantial hysteresis - the behaviour for the FCC structure 
being due to the difference in positions of T c (n) and 1/A 2 (n) maxima seen in Fig. lb. At 
present the available experimental data [ffiSSMH] correspond mainly to SQ lattices (cuprates) 
and alternating cubic lattices (Chevrel phases, doped BaBi0 3 ) and except of the overdoped 
cases they compare well with the present plots for these lattices in the intermediate coupling 
range (\U\/B < 1). 

The experimentally found deviations from the T c /T™ vs. (A m /A) 2 universal dependence 
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in the strongly overdoped regime can be explained either by the asymmetry in T>(e) (cf Fig 
9c) due to further neighbour hopping in SC and SQ structures or by the effects of intersite 
Coulomb interactions jy| , which can be different for various families of materials. 

Very recently Locquet [E3| have reported the penetration depth measurements for La2_ x 
Sr x Cu04 films as a function of doping, extending for the first time in the range from heavily 
underdoped to heavily overdoped. As we see from Fig. 10 the theoretical plots \{n) / \{n = 1) 
for our simple effective model for SQ lattice fit surprisingly well these experimental data. 
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APPENDIX A 



For rectangular DOS defined by Q47f) . one can derive from Eqs. (|10| - p!2|) analytical expressions 
determining A, H c , £ G _ L and function of n and \U\/D at T = 0. 

In particular, from Eqs.(ffl-^) taken at T = one obtains [13 

2D 



A = Jn(2-n)D/sinh 



(Al) 



fi ss = (n — 1)D coth 



2D 

W\ 



(A2) 



(n-l)D 



(A3) 



D 



U\ 



E™° = n(2-n) - 

o 2 v ' 4 



(A4) 



E, 



-n{2 - n) coth ^ - ^n 2 



D 

T 



|[/| i 

whereas 1/A 2 (0) = — K dm calculated from Eqs.(§) and 

1 A-neH 



is: 



A| ft 2 c 2 a 
Taking into account Eqs 
for if c , £ G _ L and n: 



-n(2 — n) 



2D 2D , 2£r 

coth M"M smh M, 



(A5) 



(A6) 



in Eqs.(|TTp and fli"2] ) one finds the following expressions 



# 2 



A/ 



3vr \ v a 2 (A)a ± (A) / 



(A7) 



A/ 



TTiNO 771SS 

h Q - E - 



Dn(2 - n) 



coth 



2D 

W\ 



(A8) 



SG-L 



2V2Z 



1 



coul |t/| \u\ V 1 |i/|> 



coth^ 



where D = Zt. 



(A9) 
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%V2Z 



W7T a« 



\ 



coth 2£ - 1 



ro th 25 _ 2D / • h 2D 

com \u\ \u\ V 1 |t/| 



n(2 — n) 



(A10) 



In the limits (i) \U\/D <C 1 and (ii) |t/|/-D ^> 1 the above equations take the form: 
for (i) 



1 AireH 



\ 2 r h 2 ch 



n(2 — n) 



1 



8D 



AD 



wr p \~w\ 



(All) 



Ae = Dn(2 — n)exp 



4D 



(A12) 



(A13) 



for (ii) 



$nV2Z 



7Tv 7T 



a tl sjn(2 - n) 



1 167re 2 Z t 2 



(A14) 



(A15) 



A; 



Dn(2 - n) 



\U\ 2£> ' 
2D 3|C/| 



(A16) 



G-L " 



0||t 



(A17) 



For \U\/D ^> 1 the results for if c ( and consequently also for £ G _ L and k ) are unrealistic 
as Eq° in this limit largely overestimated by MFA. As we have quoted in Sec. 2, one can 
correct these results by using Eq. ([3~8|) instead of (|A4j) in the calculation of H c . 
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FIGURE CAPTIONS 



Fig. 1 Concentration dependence of T c , the gap parameter A and 1/A 2 at T=0 K: 
\U\/B=0.1 (B= bandwidth); a) SC lattice, b) FCC lattice, c) model rectangular den- 
sity of states, (here: e 2 /h 2 c 2 a = l, a-lattice constant, 5-bandwidth) . 

Fig. 2 Concentration dependence of 1/A 2 at T=0 K for different values of U/B for SQ 
lattice. 

Fig. 3a Reduced square value of thermodynamic critical field H c at T = (h 2 c / B = 
H 2 a 2 /(8irB)) as a function of \U\/B for n— 1 for SQ lattice. HFA (o), perturbation the- 
ory: th order (obtained using Eqs. (|IT| , |l"8l , f42"|) for J = 0) (*), 2 nd order (obtained using 



Eqs.([TT||T|j4|) for (+) 



Fig. 3b h 2 c J B vs n in T = plotted for increasing values of \U\/B in HFA for SQ lattice, 
(numbers to the curves). 

Fig. 4 The G-L correlation length £ G _ L as a function of \U\/B for n = 1 for SQ lattice. HFA 
(solid line), perturbation theory: th order (obtained by putting J = in Eqs. ([T2] , p!8|j42] )) 
(dotted line), 2 nd order (obtained using Eqs.(0,|l|,||) for J ^ 0) (dashed line), T = 

Fig. 5 f G _ L (0) vs n for SQ (full circles) and FCC (triangles) lattices: \U\/B = 0.3 . 

Fig. 6 Concentration dependence of k/ko (k = A/£ G _ L at T = 0, K = K(n = l)) calculated 
for SQ lattice for several fixed values of \U\/B (numbers to the curves). The corresponding 
plot obtained for a model rectangular DOS for arbitrary \U\/B is shown by a dashed line. 

Fig. 7 The boundary between local and nonlocal electromagnetic behaviour in the \U\/B 
- n parameter space calculated from the equation k = 1/\/2 for t = 0.15 eV, Oy =3.85A and 
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a ± = 7.7A for rectangular DOS. 

Fig. 8 The values of \U\/B for the attractive Hubbard model which provide the Ginzburg 
ratio k = 10 2 . The plots as a function of n for SQ lattice, £ = 0.15 eV, ay = 3.85 A, 
a ± = 7.7A (solid line) and for rectangular DOS £ = 0.15eV, a = 7.7A (dashed line). 

Fig. 9 vs. (^) 2 at T=0 K a) for the SQ lattice, b) SC lattice ^=2 (— ), ^=0.5 (■••), 
^=0.2 (- -) and c) for FCC lattice: M = 5 (- -), ^ = 1 (- -), ^ = 0.2 (— ). A m corre- 
sponds to the maximum transition temperature T™. (o)-Tl 2 Ba2 Ca2Cu3 Oio, Tl .5 Pbo.5Sr 2 
Ca 2 Cu 3 9 , Bi 2 .xPb x Sr 2 Ca2Cu3 0i 6 ; (o)-Yi_ x Pr x Ba 2 Cu 3 6 .9 7 ; (A)-YBa 22 Cu 3 O x ; (v)-La 2 _ x 
Sr x Cu04; (*)-Bi 2 Sr 2 Cai_ x Y x Cu 2 Os+<5; (n)-LaMo6 Se 2 , PbMoeSs, SnMo6 S4Se4, SnMoeSySe, 
SnMoeSi Se7, LaMogSs, PbMog S4Se4; (©)-Tl 2 Ba 2 Cu06+5- Experimental data taken from 
22|fl- 



Fig. 10 at T=0 K for SQ lattice (^=0.1 (*), ^=0.5 (o)), compared with experi- 

mental results (full squares) |23|] La 2 _ x Sr x Cu04 films ( -39 nm) (squares), Ao = A(x= 0.159) = 
[6700±500]A. The results for a model square density of states almost coincide with the curve 
(o) for arbitrary \U/B\. 
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